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1. Introduction 

Electrical engineers and physicists are naturally very interested in the noise 
of circuits, amplifiers and detectors. This noise has many origins, some of 
which are completely unavoidable. For example, a dissipative element (a 
resistor) at finite temperature inevitably generates Johnson noise. Engi- 
neers long ago developed spectrum analyzers to measure the intensity of 
this noise. Roughly speaking, these spectrum analyzers consist of a resonant 
circuit to select a particular frequency of interest, followed by an amplifier 
and square law detector (e.g. a diode rectifier) which measures the mean 
square amplitude of the signal at that frequency. 

With the advent of very high frequency electronics operating at low 
temperatures, we have entered a new regime Hlo > k-gT, where quantum 
mechanics plays an important role and one has to begin to think about 
quantum noise and quantum-limited amplifiers and detectors. This topic is 
well-studied in the quantum optics community and is also commonplace in 
the radio astronomy community. It has recently become of importance in 
connection with quantum computation and the construction of mesoscopic 
electrical circuits which act like artificial atoms with quantized energy 
levels. It is also important for understanding the quantum measurement 
process in mesoscopic systems. 

In a classical picture, the intensity of Johnson noise from a resistor van- 
ishes linearly with temperature because thermal fluctuations of the charge 
carriers cease at zero temperature. One knows from quantum mechanics, 
however, that there are quantum fluctuations even at zero temperature, 
due to zero-point motion. Zero-point motion is a notion from quantum me- 
chanics that is frequently misunderstood. One might wonder, for example, 
whether it is physically possible to use a spectrum analyzer to detect the 
zero-point motion. The answer is quite definitely yes, if we use a quantum 



rjsdraftl2.tex; 1/02/2008; 22:49; p.l 



2 



system! Consider for example a hydrogen atom in the 2p excited state lying 
3/4 of a Rydberg above the Is ground state. We know that this state is 
unstable and has a lifetime of only about 1 ns before it decays to the ground 
state and emits an ultraviolet photon. This spontaneous decay is a natural 
consequence of the zero-point motion of the electromagnetic fields in the 
vacuum surrounding the atom. In fact, the rate of spontaneous decay gives 
a simple way in which to measure this zero point motion of the vacuum. 
Placing the atom in a resonant cavity can modify the strength of the noise 
at the transition frequency, and this effect can be measured via a change 
in the decay rate. 

At finite temperature, the vacuum will contain blackbody photons which 
will increase the rate of decay due to stimulated emission and also cause 
transitions in the reverse direction, Is — > 2p, by photon absorption. With 
these ideas in mind, it is now possible to see how to build a quantum 
spectrum analyzer. 

The remainder of this article is organized as follows. First we describe 
the general concept of a two-level system as a quantum spectrum analyzer. 
We next review the Caldeira-Leggett formalism for the modelling of a 
dissipative circuit element, such as a resistor, and its associated quantum 
noise. Then, a brief discussion of the single Cooper-pair box, a circuit which 
behaves as a two- level system or qubit, is given. We then discuss the effects 
of a dissipative electromagnetic environment on the box, and treat the case 
of a simple linear, but nonequilibrium environment, consisting of a classical 
tunnel junction which produces shot noise under bias. Finally, we describe 
a theoretical technique for calculating the properties of a Cooper-pair box 
coupled to a measurement system, which will be a nonlinear, nonequilibrium 
device, such as a single-electron transistor. Equivalently, this allows one 
to calculate the full quantum noise spectrum of the measurement device. 
Results of this calculation for the case of a normal SET are presented. 



2. Two- level systems as spectrum analyzers 

Consider a quantum system (atom or electrical circuit) which has its two 
lowest energy levels eo and t\ separated by energy Eq\ = Huioi- We suppose 
for simplicity that all the other levels are far away in energy and can be 
ignored. The states of any two-level system can be mapped onto the states 
of a fictitious spin- 1/2 particle since such a spin also has only two states 
in its Hilbert space. With spin up representing the ground state (\g)) and 
spin down representing the excited state (|e)), the Hamiltonian is (taking 
the zero of energy to be the center of gravity of the two levels) 

rj ^01 m 
H = 7T- °z- (1) 
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In keeping with the discussion above, our goal is to see how the rate of 'spin- 
flip' transitions induced by an external noise source can be used to analyze 
the spectrum of that noise. Suppose for example that there is a noise source 
with amplitude f(t) which can cause transitions via the perturbation 1 

V = Af(t)a x , (2) 

where A is a coupling constant. The variable f(t) represents the noise 
source. We can temporarily pretend that / is a classical variable, although 
its quantum operator properties will be forced upon us very soon. For now, 
only our two-level spectrum analyzer will be treated quantum mechanically. 

We assume that the coupling A is under our control and can be made 
small enough that the noise can be treated in lowest order perturbation 
theory. We take the state of the two-level system to be 

i*w>-(2w)- (3) 

In the interaction representation, first-order time-dependent perturbation 
theory gives 

\Mt)) = \m)-y*dTV { T)\m)- (4) 

If we initially prepare the two-level system in its ground state then the 
amplitude to find it in the excited state at time t is 

i A f* 

a e = -— J dr (e\a x (r)\g)f(r) + 0(A 2 ), (5) 

f dr e^ 01T f{r) + O(A 2 ). (6) 
Jo 

We can now compute the probability 

p e (t) = |a e | 2 = ^J*J*dT 1 dr 2 e- iU01 ^-^f(r 1 )f(r i ) + O(A 3 ) (7) 
We are actually only interested on the average time evolution of the system 



iA '"' 



fr 2 Jo Jo 



Pe(t) = '- W l I dndne-^ ^^)(/(r 1 )/(r 2 )) + 0(A 3 ) (8) 



1 The most general perturbation would also couple to a y but we assume that (as is 
often, though not always, the case) a spin coordinate system can be chosen so that the 
perturbation only couples to a x . Noise coupled to a z commutes with the Hamiltonian 
but is nevertheless important in dephasing coherent superpositions of the two states. We 
will not discuss such processes here. 
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We can now perform a change of variables in the integrals, r = t\ — T2 and 
T = (ti + T2) /2, and we get 

42 /•* /-B(T) 

p c (t) =-t dT dr (/(T + r/2)/(T _ r/2)) + 0(i4 3 } (g) 

ft JO J-B(T) 

where 

B(T) = T if T < t/2 

= t - T if T > t/2. 

Let us now suppose that the noise correlation function is stationary (time 
translation invariant) and has a finite but small autocorrelation time Tf. 
Then for t S> 17 we can set the bound B (T) to infinity in the last integral 
and write 

A2 r-t noo 

p e (t) = —^dTjjr </(r)/(0)> + 0(A 3 ) (10) 

The integral over r is effectively a sum of a very large number N ~ i/ry of 
random terms 2 and hence the value undergoes a random walk as a function 
of time. Introducing the noise spectral density 

/+00 
dr e ^(/(r)/(0)>, (11) 
-00 

we find that the probability to be in the excited state increases linearly 
with time, 3 

Pe(t) = t^Sf(-Uoi) (12) 

The time derivative of the probability gives the transition rate 

r T = ^Sfi-uu) (13) 

Note that we are taking in this last expression the spectral density on the 
negative frequency side. If / were a strictly classical source (/(t)/(0)) would 



2 The size of these random terms depends on the variance of / and on the value of 
u)oiTf For LuoiTf S> 1 the size will be strongly reduced by the rapid phase oscillations of 
the exponential in the integrand. 

3 Note that for very long times, where there is a significant depletion of the probability 
of being in the initial state, first-order perturbation theory becomes invalid. However, for 
sufficiently small A, there is a wide range of times r/ < ( « 1/r for which Eq. 12 is valid. 
Eqs. 13 and 14 then yield well-defined rates which can be used in a master equation to 
describe the full dynamics including long times. 
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be real and Sf(— woi) = Sf(+u>oi). However, because as we discuss below 
/ is actually an operator acting on the environmental degrees of freedom, 
[/(r), /(0)] + and S f (-io i) + S f (+u 01 ). 

Another possible experiment is to prepare the two-level system in its 
excited state and look at the rate of decay into the ground state. The 
algebra is identical to that above except that the sign of the frequency is 
reversed: 

r i = ^(+^01). (14) 

We now see that our two-level system does indeed act as a quantum spec- 
trum analyzer for the noise. Operationally, we prepare the system either 
in its ground state or in its excited state, weakly couple it to the noise 
source, and after an appropriate interval of time (satisfying the above in- 
equalities) simply measure whether the system is now in its excited state or 
ground state. Repeating this protocol over and over again, we can find the 
probability of making a transition, and thereby infer the rate and hence the 
noise spectral density at positive and negative frequencies. Note that in con- 
trast with a classical spectrum analyzer, we can separate the noise spectral 
density at positive and negative frequencies from each other since we can 
separately measure the downward and upward transition rates. Negative 
frequency noise transfers energy from the noise source to the spectrometer. 
That is, it represents energy emitted by the noise source. Positive frequency 
noise transfers energy from the spectrometer to the noise source. 4 In order to 
exhibit frequency resolution, Auj, adequate to distinguish these two cases, 
it is crucial that the two-level quantum spectrometer have sufficient phase 
coherence so that the linewidth of the transitions satisfies the condition 
luoi/ Alu > maxffceT'/^woi) !]• 

In thermodynamic equilibrium, the transition rates must obey detailed 
balance T^/T^ = e^^ 01 in order to give the correct equilibrium occupancies 
of the two states of the spectrometer. This implies that the spectral densities 
obey 

S f (+uoi) = e^S f {-u Jm ). (15) 

Without the crucial distinction between positive and negative frequencies, 
and the resulting difference in rates, one always finds that our two level 

4 Unfortunately, there are several conventions in existence for describing the noise 
spectral density. It is common in engineering contexts to use the phrase 'spectral density' 
to mean Sf(+ui) + Sf(-iv). This is convenient in classical problems where the two 
are equal. In quantum contexts, one sometimes sees the asymmetric part of the noise 
Sf(+uj) — Sf(-uj) referred to as the 'quantum noise.' We feel it is simpler and clearer to 
simply discuss the spectral density for positive and negative frequencies separately, since 
they have simple physical interpretations and directly relate to measurable quantities. 
This convention is especially useful in non-equilibrium situations where there is no simple 
relation between the spectral densities at positive and negative frequencies. 
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systems are completely unpolarized. If, however, the noise source is an 
amplifier or detector biased to be out of equilibrium, no general relation 
holds. 

We now rigorously treat the quantity /(t) as quantum operator in the 
Hilbert space of the noise source. The previous derivation is unchanged, and 
Eqs. (13,14) are still valid provided that we interpret the angular brackets 
in Eq. (8) as representing the quantum statistical expectation value for the 
operator correlation (in the absence of the coupling to the spectrometer) 

/+oo 
dre^^p QQ (a|/(r)| 7 )( 7 |/(0)|a) (16) 

where for simplicity we have assumed that (in the absence of the coupling 
to the spectrometer) the density matrix is diagonal in the energy eigen- 
basis and time-independent (but not necessarily given by the equilibrium 
expression). This yields the standard quantum mechanical expression for 
the spectral density 

/+oo 
dTe^^Paae*^ >*|M/l7>| 2 (17) 

= 2^hY^Pa a \{a\f\l)\ 2 5{e^-e a -nw). (18) 

a, 7 

Substitution of this into Eqs. (13,14) we derive the familiar Fermi Golden 
Rule expressions for the two transition rates. 

In standard courses, one is not normally taught that the transition 
rate of a discrete state into a continuum as described by Fermi's Golden 
Rule can (and indeed should!) be viewed as resulting from the continuum 
acting as a quantum noise source. The above derivation hopefully provides 
a motivation for this interpretation. 

One standard model for the continuum is an infinite collection of har- 
monic oscillators. The electromagnetic continuum in the hydrogen atom 
case mentioned above is a prototypical example. The vacuum electric field 
noise coupling to the hydrogen atom has an extremely short autocorre- 
lation time because the range of mode frequencies to a (over which the 
dipole matrix element coupling the atom to the mode electric field E a 
is significant) is extremely large, ranging from many times smaller than the 
transition frequency to many times larger. Thus the autocorrelation time of 
the vacuum electric field noise is considerably less than 10 _15 s, whereas the 
decay time of the hydrogen 2p state is about 10~ 9 s. Hence the inequalities 
needed for the validity of our expressions are very easily satisfied. 

Of course in the final expression for the transition rate, energy con- 
servation means that only the spectral density at the transition frequency 
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enters. However, in order for the expression to be valid (and in order for 
the transition rate to be time independent), it is essential that there be a 
wide range of available photon frequencies so that the vacuum noise has an 
autocorrelation time much shorter than the inverse of the transition rate. 

3. Quantum Noise from a Resistor 

Instead of an atom in free space, we might consider a quantum bit capaci- 
tively coupled to a transmission line. The transmission line is characterized 
by an inductance per unit length I and capacitance per unit length c. A 
semi-infinite transmission line presents a frequency-independent impedance 
Z = Rq = \j 1 tjc at its end and hence acts like an ideal resistor. The 
dissipation is caused by the fact that currents injected at one end launch 
waves which travel off to infinity and do not return. Very conveniently, 
however, the system is simply a large collection of harmonic oscillators (the 
normal modes) and hence can be readily quantized. This representation of 
a physical resistor is essentially the one used by Caldeira and Leggett [1] 
in their seminal studies of the effects of dissipation on tunneling. The only 
difference between this model and the vacuum fluctuations in free space 
discussed above is that the relativistic bosons travel in one dimension and 
do not carry a polarization label. This changes the density of states as a 
function of frequency, but has no other essential effect. 
The Lagrangian for the system is 



where j(x, t) is the local current density and q(x,t) is the local charge den- 
sity. Charge conservation connects these two quantities via the constraint 



in terms of which the current density is j(x, t) = —dfO(x, t) and the charge 
density is q(x,t) = d x 6(x,t). For any well-behaved function 0(x,t), the 
continuity equation is automatically satisfied so there are no dynamical 
constraints on the 6 field. In terms of this field the Lagrangian becomes 




(19) 



d x j(x,t) + d t q(x,t) = 0. 



(20) 



We can solve this constraint by defining a new variable 




(21) 




(22) 
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The Euler-Lagrange equation for this Lagrangian is simply the wave equa- 
tion v 2 d^9 — d 2 = where the mode velocity is v = 1/VIc. 

^From Eq. (21) we can deduce that the proper boundary conditions 
(in the absence of any coupling to the qubit) for the 9 field are 6(0, t) = 
0(L, t) = 0. (We have temporarily made the transmission line have a fi- 
nite length L.) The normal mode expansion that satisfies these boundary 
conditions is 



^*) = \/l&«W sill X' (23) 

n=l 

where (p n is the normal coordinate and k n = Substitution of this form 
into the Lagrangian and carrying out the spatial integration yields a set of 
independent harmonic oscillators representing the normal modes. 

£ = E^n) 2 -±- c k^l (24) 
n=l 

From this we can find the momentum p n canonically conjugate to (p n and 
quantize the system to obtain an expression for the voltage at the end of the 
transmission line in terms of the mode creation and destruction operators 

v-yflw<M)-ii>^(4+*.). (*) 

n=l 

The spectral density of voltage fluctuations is then found to be 
2 °° Ml 

S v (lj) =2TT-^2^{n 1 {hfl n )5(uj + n n ) + [n 1 (hn n ) + l}5(uj-n n )}, (26) 

n=l 

where n^{hJ) is the Bose occupancy factor for a photon with energy hio. 
Taking the limit L — > oo and converting the summation to an integral yields 

S v {u) = 2R h\u\{n J {H\uj\)e(-u) + [n^hio) + l]0(w)}, (27) 

where is the step function. We see immediately that at zero temperature 
there is no noise at negative frequencies because energy can not be extracted 
from zero-point motion. However there remains noise at positive frequencies 
indicating that the vacuum is capable of absorbing energy from the qubit. 

A more compact expression for this 'two-sided' spectral density of a 
resistor is 

S "M = T^kf < 28 > 

which reduces to the more familiar expressions in various limits. For ex- 
ample, in the classical limit ksT S> Mo the spectral density is equal to the 
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Johnson noise result 5 

S v {oj) = 2R k B T, (29) 
which is frequency independent, and in the quantum limit it reduces to 

S v (u) = 2R nwe(uj). (30) 

Again, the step function tells us that the resistor can only absorb energy, 
not emit it, at zero temperature. 

If we use the engineering convention and add the noise at positive and 
negative frequencies we obtain 

S v (u) + S v (-u) = 2R hvcoth-J^- (31) 

for the symmetric part of the noise, which appears in the quantum fluctuation- 
dissipation theorem [2]. The antisymmetric part of the noise is simply 

S v (u) - Sv(-u) = 2R friu. (32) 

This quantum treatment can also be applied to any arbitrary dissipative 
network [3]. If we have a more complex circuit containing capacitors and 
inductors, then in all of the above expressions, i?o should be replaced by 
ReZ(uj) where Z(u) is the complex impedance presented to the qubit. 



4. The Single Cooper-Pair Box: a Two-Level Quantum Circuit 

The Cooper-pair box (CPB) is a simple circuit [4], consisting of a small 
superconducting "island" , connected to a large reservoir via a single small- 
capacitance Josephson junction, depicted as a box with a cross (Fig. 1). The 
island is charge biased by applying a voltage (V g ) to a nearby lead, called 
the gate, which has a small capacitance to the island, C g . The junction 
is characterized by its capacitance, Cj, and its tunnel resistance, Rj. At 
temperatures well below the transition temperature of the superconductor 
(Tq ~ 1.5 K for the usual Al/AlOx/Al junctions), none of the many (~ 10 9 ) 
quasiparticle states on the island should be thermally occupied, and the 
number of Cooper-pairs on the island is the only relevant degree of freedom. 

We may then write the Hamiltonian for the box in terms of the states of 
different numbers of pairs on the island, which are eigenstates of the number 
operator, h\n) = n\n). The box Hamiltonian consists of an electrostatic 
term, plus a Josephson term describing the coupling of the island to the 

5 Note again that in the engineering convention this would be Sv(u) = 4RoksT. 
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Figure 1. a) Circuit diagram of Cooper-pair box. b) Pseudo-spin representation of the 
energies of Cooper-pair box. The density matrix for the two pure eigenstates lie along 
the total effective field, collinear with the z' axis. 



lead, 



" — H electrostatic ~\~ H j osephs 



AE C X> 



\n){n\ 



Ej 



^2(\n + l)(n\ + h.c.) 



(33) 
(34) 



The energy scale for the electrostatic interaction is given by the charging 
energy, Eq = e 2 /2Cs, where Cs = Cj + C g is the total island capacitance, 
while the Josephson energy, Ej, is set by the tunnel resistance and the gap 
of the superconductor, 



Ej 



hA 

8e 2 Ri 



ARk 

8 Ra 



(35) 



The electrostatic term is easily modulated by changing the voltage on the 
gate; the quantity n g = C g V g /2e that appears in the Hamiltonian corre- 
sponds to the total polarization charge (in units of Cooper pairs) injected 
into the island by the voltage source. 

This Hamiltonian leads to particularly simple behavior in the charge 
regime, when the electrostatic energy dominates over the Josephson cou- 
pling, 4Ec 3> Ej. In this case we can restrict the discussion to only two 
charge states, |n = 0) and |n = 1). For convenience we can reference the 



energies of the two states to their midpoint, E m id 
the Hamiltonian now becomes 



H 



-E e i —Ej 
-Ej E e l 



4E c (l-2n g y, so that 



(36) 



where E e i is the electrostatic energy that is now linear in the gate charge, 
E e i = ^Ec(l—2n g ). It is also now apparent that the Hamiltonian is identical 
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Figure 2. Energies (left) of ground and excited states of a Cooper-pair box with 
Ec — Ej vs. dimensionless gate charge, n g = C g V g /2e. The expectation value of a 
charge measurement, (n), (right) for the ground (solid line) and excited (dotted line) 
states vs. n g . 



to that of a fictitious spin-1/2 particle, 

tj E el Ej 

H = — 2~<?z ~ yffx, (37) 

under the influence of two psuedo-magnetic fields, B z = E e i and B x = Ej, 
as depicted in Fig. 1. In other words, the box is a qubit or two-level system 6 . 
The state of the system is in general a linear combination of the states 
\n = 0} and |n = 1). The state can be depicted using the density matrix, 
which corresponds to a point on the Bloch sphere, where the north pole 
(+z-direction) corresponds to |n = 0). The ground and excited states of 
the system will be aligned and anti-aligned with the total fictitious field, 
i.e. in the ±z' directions. 

It is also apparent from this discussion that the states of the box can 
be easily manipulated by changing the gate voltage. The energies of the 
ground and excited states, as a function of n g , are displayed in Figure 2. 
The energy difference between the ground and excited bands varies from 
AEc at n g = 0, 1, to a minimum at the charge degeneracy point, n g = 1/2. 
At this point, the Josephson coupling leads to an avoided crossing, and the 
splitting is Ej. 

Also plotted is the expectation value of the number operator, (|n|), 
which is proportional to the total charge on the island. In the geometrical 

6 Of course, this is an approximation, as there are other charge states (|n = 2), etc.) 
which are possible, but require much higher energy. Even outside the charge regime (i.e. 
when Ej > AEc) the two lowest levels of the box can be used to realize a qubit [5]. In 
this case, the two states do not exactly correspond with eigenstates of charge, and matrix 
elements are more complicated to calculate. Nonetheless, this regime can also be used as 
an electrical quantum spectrum analyzer. 
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picture of Fig. lb, a measurement of charge (n) is equivalent to projecting 
the state on the z-axis, h = ^(1 — a z ). We see that as the gate charge 
is changed from to 1, the ground state is initially |n = 0), and the 
character of the ground and excited states interchange on passing through 
the degeneracy point, leading to the transition between (n) = and 1, which 
is broadened by quantum fluctuations (the a x coupling) . At the degeneracy 
point, the ground and excited states lie in the ±x-directions, i.e. they are 
symmetric and antisymmetric combinations of the two charge states. In 
general, we will denote the ground and excited state of the CPB at a 
particular gate voltage as \g) and |e), which are given in terms of the charge 
states by \g) = cos(0/2)|O)+sin(0/2)|l) and|e) = - sin(0/2)|O)+cos(0/2)|l) 
respectively, where = &ictan\Ej / E e i] is a function of the gate voltage. 

A nice property of the CPB in this regime is that the various matrix 
elements can be calculated in a straightforward way. For example, the 
expectation value of n in the ground state, {g\h\g), is therefore equal to 
1/2(1 — (g\<7 z \g)) = sin 2 (0/2), from which we can find the ground state 
charge as shown in Fig. 2. A perturbation in the gate charge, due for 
example to a fluctuation or change in the applied gate voltage, will lead 
to a proportional change in the electrostatic energy, or the z-component of 
the fictitious magnetic field. Such a perturbation will cause both dephasing 
and transitions between states. 

5. General Discussion of CPB Coupled to a Dissipative Environment 

In the previous section we described the Hamiltonian and the eigenstates for 
a Cooper-pair box which is "charge-biased," i.e. controlled with a voltage 
applied to a gate capacitor C g , as shown in Fig. 1. In our earlier treatment 
of the box, the voltage and the dimensionless gate charge, n g were treated 
as fixed parameters of the Hamiltonian (c-numbers). In this case, the box's 
evolution is purely deterministic and conservative. However, it is impossible, 
even in principle, to control such a voltage with arbitrary precision at all 
frequencies. In Fig. 1, the idealized source of the gate voltage is drawn in 
series with an impedance Z{oj) of the gate lead. Generally, this gate lead 
will be connected to external wiring (a transmission line), with a typical 
real impedance comparable to the impedance of free space (~ 50 VI) at 
the microwave transition frequencies of the box. ^From the fluctuation- 
dissipation theorem we know that this impedance will introduce noise on 
the gate voltage, even at zero temperature. 

There are several effects of the voltage noise on the box, or the coupling 
of our spin- 1/2 circuit to the many external degrees of freedom represented 
by the gate impedance. First, even at zero temperature, we will find a finite 
excited state lifetime, Tl, for the box. Second, at finite temperature, we will 
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find a finite polarization of our psuedo-spin, i.e. some steady-state proba- 
bility to find the spin in its excited state. Finally, the gate noise introduces 
a random effective field felt by the spin, and a loss of phase coherence for a 
superpostition state. It is this last effect which is most important in making 
high-fidelity qubits and performing quantum computations, but it is the 
first two which depend most explicitly on the quantum nature of the noise. 
We deal in this manuscript with only these first two features of the box's 
coupling to the electromagnetic environment, and ignore the dephasing 7 . Of 
course, the other parameter in the Hamiltonian, the Josephson energy, can 
in principle fluctuate, especially as in many experiments the box's junction 
is split into a small SQUID in order to provide external tuning of Ej with 
an applied flux. We concentrate here only on the voltage noise (fluctuations 
in the a z part of the Hamiltonian) for simplicity. 

We begin with a very simple treatment of the dynamics of the two-level 
system under the influence of the gate voltage noise. We are interested in 
the ensemble-averaged behavior of our psuedo-spin, which is best described 
using the density matrix approach, and is detailed in Section 7 on the 
coupling of the box to a measuring SET. The basic effects, apart from 
dephasing, however, can be captured simply by examining the probabilities 
p g and p e of finding the box in its ground (\g)) or excited (|e)) states. 
The noise of the external environment can drive transitions from ground to 
excited state and vice- versa, at rates and T^, respectively. The coupled 
master equations for these probabilities are 



Of course conservation of probability tells us that p e + p g = 1, so we 
introduce the polarization of the spin-1/2 system, P = p g — p e . In steady- 
state, the detailed balance condition is p e T^ = p g Tp The two rates T| 
and 1^ are related by Equations 13 and 14 to the spectral densities of the 
noise at negative and positive frequencies. We see immediately that if the 
spectral density is symmetric (classical!), then the rates for transitions up 
and down are equal, the occupancies of the two states are exactly equal, and 
the polarization of the psuedo-spin is identically zero. It is the quantum, or 
antisymmetric, part of the noise which gives the finite polarization of the 
spin. Even in NMR, where the temperature is large compared to the level 
splitting (Itlooi < k^T), this effect is well-known and crucial, as the small 
but non-zero polarization is the subject of the field! 



7 For a nice recent treatment of dephasing in Josephson junctions, see Ref. [6]. 



dpe 

dt 

dpg_ 

dt 



PgT^ - Pe Y i 



(38) 



(39) 
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Solving for the steady-state polarization, we find 



£| - T T = S(+u 01 ) - S(-uj 01 ) 
r^ + Tt ~ 5(+w i) + S'(-woi) 



(40) 



An measurement of the steady-state polarization allows one to observe the 
amount of asymmetry in the noise, so the two-level system is a quantum 
spectrum analyzer. 

If we can create a non-equilibrium polarization, P = P ss + AP (a pure 
state is not necessary) of our two-level system, we expect it to return to the 
steady state value. Substituting the modified probabilities p e (t) = p ess — 
AP(t)/2 and p g (t) = p 9ss + AP(i)/2 into our master equations above, we 
find an equation for the deviation of the polarization 



Thus the system decays to its steady-state polarization with the relaxation 
rate Ti = I\ + r\ = (A/h) 2 [S(— loqi) + S(+ujqi)} related to the total 
noise at both positive and negative frequencies. In NMR, the time 1/Ti 
is referred to as XI. In the zero-temperature limit, there is no possibility of 
the qubit absorbing energy from the environment, so = 0, and we find 
full polarization P = 1, and a decay of any excited state population at a 
rate = 1/TI which is the spontaneous emission rate. 

It is worth emphasizing that a quantum noise source is always char- 
acterized by two numbers (at any frequency), related to the positive and 
negative frequency spectral densities, or to the symmetric and antisym- 
metric parts of the noise. These two quantities have different effects on a 
two-level system, introducing a finite polarization and finite excited-state 
lifetime. Consequently, a measurement of both the polarization and Tl of a 
two-level system is needed to fully characterize the quantum noise coupled 
to the qubit. Such measurements in electrical systems are now possible, 
and some of us [7] have recently performed such a characterization for 
the specific case of a CPB coupled to a superconducting single-electron 
transistor. 

Our discussion in this section uses the language of NMR to describe 
the effects on the two-level system. There are, however, several possible 
protocols 8 for measuring the quantum noise, and several different "basis 
sets" or measured quantities which describe the noise process or the quan- 
tum reservoir to which the two-level system is coupled. Table 1 contains 
a "translation" between the specific pairs of quantities that are commonly 

8 The idea of watching the decay from the pure states |e) and \g) to measure 
SV(±|k>oi|) separately was described in Section 2. 



d(AP(t)) 



= -AP(t)(r T + r i ). 



(41) 



dt 
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TABLE I. Different ways to characterize a quantum reservoir. 



Fermi Golden Rule r T (w) = ^Sv(+\uj\) I^w) = j?Sv(-\uj\) 



Fluct.-Diss. Relation n 7 (w) = 2r T /(r| - r T ) Re[Z(u)] = ^(Ij - T T ) 



NMR 



ri = (r i + r T ) 



-i 



^ = (ri-r T )/(r i+ r T ) 



Quantum Optics 




Aev 



'instein 



used in different disciplines, and their relation to the positive and negative 
noise spectral densities. In all cases, though, two separate numbers are 
required to specify the properties of a quantum reservoir. 

6. The Box Coupled to an Ohmic Environment 

We can now proceed to the effects of a specific dissipative coupling to the 
Cooper-pair box. The noise on the gate voltage will lead to a fluctuation of 
the gate charge parameter, n g , and thus to a fluctuation of the electrostatic 
energy, i.e. the a z term in the Hamiltonian (Eq. 37). Depending on the 
average value of n g , this fluctuation will consist of fluctuations which arc 
both longitudinal (|| to a' z ) and transverse (_L to a' z ). To calculate the rates 
of transitions between the states |e) and \g), we need to find the coupling 
strength A of this perturbation in the a x direction. Referring to Fig. 1, 
we see that a z = cos(6)a' z — sin(#)cr^.. If we let the gate charge now be 
n g (t) = fig + 5n g (t), we may rewrite the Hamiltonian Eq. 37 in the new 
eigenbasis as 



H = — z^V + 4E C cos(e)5n g (t)a' z - 4E C sin(0) 5n g {t)a' x . (42) 



The time-varying term in the o~' z direction will effectively modulate the 
transition frequency, loqi = Eqi/H, and cause dephasing. In terms of the 
gate voltage noise, V(t), the a' x perturbation term has the form AV(t)a' x = 
eKsm(9)V(t)a' x , where e is the electron's charge and k = C 9 /Cy, is the 
capacitive coupling. Using Eq. 14, we find 



If the environment is effectively at zero temperature (fajoi ^ &bT), then 
Sv{+wqi) = 2huiQiR, and the quality factor of the transition is 




(43) 



Q = woi/Tj, 



_L Rk_ 

? sin 2 9 47ri? 



(44) 
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where Rk = h/e 2 is the resistance quantum. 

For a finite temperature, we have rates in both directions, and the 
polarization of the psuedo-spin is given by the ratio of the antisymmetric 
(Eq. 32) to symmetric (Eq. 31) spectral densities, P = tanh(fia;oi/2fcB7 1 ), 
as one expects for any two-level system at temperature T. An example of 
the average charge state of a Cooper-pair box at finite temperature, and of 
the polarization and equilibration time Tl, are shown in Figure 3. As the 
gate voltage is varied, the transition frequency of the box changes from a 
maximum near n g = 0, to a minimum c^oi = Ej jh at the degeneracy point 
n g = 0.5 and then back again. We see that the states of the box are generally 
most "fragile" near the avoided crossing. First, the energy splitting is a min- 
imum here, leading to the lowest polarization of the psuedo-spin. Second, 
because the eigenstates \g) and |e) point in the a x direction, the matrix 
elements for the voltage fluctuations of the environment are maximal, i.e. 
the noise is orthogonal to the spin. This also implies that the dephasing 
effects are minimal at this degeneracy point, which offers great advantages 
for improving the decoherence times [5], but the excited state lifetimes are 
smallest at this point. One also sees that the lifetimes become large away 
from the degeneracy, where the electrostatic energy dominates over the 
Josephson energy, which offers advantages when measuring the charge state. 
In the limit that Ej could be suppressed to zero during a measurement, the 
matrix elements (for voltage noise) vanish, and a quantum non-demolition 
(QND) measurement [8] could be performed. The idea of using the qubit 
as a quantum spectrum analyzer is precisely the reverse, where we measure 
the "destruction" in the two-level system (i.e. inelastic transitions caused 
by the coupling of the states to the environment), in order to learn about 
the quantum noise spectrum of the fluctuations. 

The Cooper-pair box can of course also be used to measure the more in- 
teresting spectral densities of nonequilibrium devices. The simplest example 
is to replace the gate resistance by a tunnel junction. If we arrange to bias 
the junction using, e.g. an inductor, a dc current / and an average dc voltage 
V can be maintained across the junction. Classically, the current noise of 
such a tunneling process is frequency independent, Si = 2el. The voltage 
noise density presented to the CPB's gate would then be Sy = 2eIR^, 
where Rt is the junction tunnel resistance. In fact, this "white" spec- 
tral density can only extend up to frequencies of order uj = eV/h, the 
maximum energy of electrons tunneling through the junction. The correct 
high-frequency form of the symmetrized noise density was calculated by 
Rogovin and Scalapino [9], 



S v H = R(hw + eV) coth 



hw + eV 



2k B T 



R (hcu - eV) coth 



hijj — eV 
2k B T 



(45) 
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T1 (jis) 



Figure 3. The box coupled to an equilibrium, Ohmic environment, i.e. a resistor, a) 
Two-sided noise spectral density, of the voltage, SV(w), for a resistor at T=0 (solid 
line) and finite temperature (dashed line) b) Average charge of box with Ec = IK, 
E.j = 0.5K when coupled with strength k = O.Of to a resistor with resistance R — 50 Q 
and temperature T = 0.5K. c) Polarization (dotted line) and relaxation time Tl for the 
same parameters. Full line is the rate of spontaneous emission, i.e. Tl at zero temperature. 



and was indirectly measured in a mesoscopic conductor using a conventional 
spectrum analyzer [10]. This noise can also be expressed [11] in its two-sided 
form 

{hw + eV)R T (hio-eV)R T 

<VK> = hui+eV H hu ,_ eV , (46) 

1 — e k B T 1 — e k B T 

and is displayed in Fig. 4. Notice that the antisymmetric part of this noise is 
the same as that of the ordinary resistor, Sv(+w) — Sv(-uj) = 2TvjjRt, and 
is independent of the voltage. Also shown in Fig. 5 is the polarization and 
relaxation time, Tl, of a CPB coupled to a shot noise environment. We see 
that full polarization is achieved only when ftwrji S> eV. For low transition 
frequencies, the polarization is inversely proportional to the current through 
the junction. Aguado and Kouwenhoven [11] have described the use of a 
double quantum dot as a two-level system to probe this behavior of the 
shot noise. 

Given our discussion so far, it is now interesting to ask what the effects 
of a real quantum measurement on a quantum circuit will be. A quantum 
measurement device will in general be neither linear, Ohmic, nor equilib- 
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Figure 4- The box coupled to an nonequilibrium, Ohmic environment, i.e. a 50 Q 
tunnel junction, a) Two-sided noise spectral density, of the voltage, «SV(w), for a junction 
at T=0.02 K, with zero voltage (dotted line), and increasing bias voltages (solid and 
dashed lines), b) Average charge of box with Ec = IK, Ej — 0.5K when coupled with 
strength n = 0.01 to a junction biased at eV = 1.5 K. c) Polarization (dotted line) and 
relaxation time Tl for the same parameters (solid) and for T=0, V=0 (dashed line). 



rium. Obviously, if we hope to characterize this measurement process, and 
to understand what one will observed when the qubit is coupled to the 
noise processes of the measuring device, we will need to calculate the full 
quantum (two-sided!) noise spectrum of the amplifier or detector. 

7. Single-Electron Transistor Coupled to a Two-Level System 

We have seen in previous sections that a two-level system (TLS) may be 
used as a "spectrum analyzer" to measure quantum noise. Here, we use this 
technique to theoretically calculate quantum noise. Instead of simply study- 
ing the "noisy" system of interest in isolation, one can study a composite 
system consisting of the "noisy" system coupled to a TLS; by calculating 
the relaxation and excitation rates of the TLS, one can efficiently calcu- 
late the quantum noise of interest 9 . We demonstrate the usefulness of this 

9 Note that the spirit of our approach is similar to that employed in the theory of full 
counting statistics [12]. There too one attaches an auxiliary spin 1/2 to the scattering 
system of interest, and studies the dynamics of the fully coupled system to obtain the 
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technique by outlining a calculation of the quantum charge noise of a single 
electron transistor (SET). This is an important example, as when an SET is 
used as an electrometer, it is this noise which determines the measurement 
backaction. 

The SET consists of a metallic island attached via tunnel junctions to 
source and drain reservoirs. It is described by the Hamiltonian: 

Hset = ( £ k-Va)c{ a c ka + E c (n-N) 2 + H T (47) 

k,a=L,R,I 

H T = t (^Iv + Zi.c.) (48) 

k,q, a=L,R 

The first term in Hset describes the kinetic energy of electrons in the leads 
(a = L,R) and on the island (a = I). The second term is the Coulomb 
charging energy which depends on n, the number of excess electrons on 
the island. This interaction term can be tuned by changing the voltage 
on a nearby gate electrode which is capacitively coupled to the island; N 
represents the dimensionless value of this voltage. Finally, Ht describes the 
tunneling of electrons through the two SET tunnel junctions; the conduc- 
tance of each junction (in units of e 2 /h) is given by g = 4ir 2 t 2 i/Q, with uq 
being the density of states. is an auxiliary operator which increases n 
by one: [n, F^ = F^ . For simplicity, we assume that the two junctions of 
the SET are completely symmetric (i.e. equal junction capacitances and 
dimensionless conductances). 

Throughout this section, we will be interested in the regime of sequential 
tunneling in the SET, where transport involves energy-conserving transi- 
tions between two charge states of the SET island, say n = and n = 1. 
These transitions are described by simple rates, which can be derived via 
Fermi's Golden rule: 

rS± 1>B = 7([AiC ± i,J (49) 

AE% ±lin = T2E C (n±±-M)±(±- 5 a , R ^j eV DS (51) 

r" ±1 n is the tunneling rate from the charge state n to n ± 1 through 
junction a; AE is the energy gained in making the tunneling transition, 
and includes contributions both from the drain-source voltage Vds an d 
from the charging energy. Sequential tunneling is the dominant transport 
mechanism when the junction conductances are small (i.e. g/(2ir) <C 1), 



statistics of charge transfer in the scatterer. 
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and the dimensionless gate voltage N is not too far away from a charge 
degeneracy point. Sequential tunneling is the most important regime for 
measurement applications, as it yields the largest SET currents. 

At low temperatures, only tunnel events which follow the voltage are 
possible. There are thus there only two relevant rates: n = — > 1 transitions 
occur through the left junction at a rate Ff , while n = 1 — > occur through 
the right junction at a rate T^. The average current will be given by: 

(I) = e^§n (52) 

10 01 

We are interested in calculating Sq(lo), the quantum noise associated 
with fluctuations of the charge on the central island of the SET. It is defined 
as: 

/oo 
rft(n(*)n(0))e- <hrt (53) 
-oo 

Note that we can equivalently think of Sq as describing the voltage noise 
of the island, as Vi s i an d = en/C^;, where Cs is the total capacitance of the 
island. In two limits, the form of the island charge noise can be anticipated. 
For lo — ► 0, the noise will correspond to classical telegraph noise- the island 
charge n fluctuates between the values and 1, with Poisson-distributed 
waiting times determined by the rates Tf and T^. We thus expect a 
symmetric, Lorentzian form [13] for the noise at low frequencies: 

For large frequencies \uj\ 3> Ec, we expect that correlations due to the 
charging energy will have no influence on the noise. The system will effec- 
tively look like two tunnel junctions in parallel, and we can use the results 
of Sec. 6 for the corresponding voltage noise. Noting that each junction 
effectively consists of a resistor and capacitor in parallel, we have at zero 
temperature: 

C 2 C 2 
S q (lo) = -f x S v (u>) -► [2hwRe Z tot Q{u)] (|w| » E c ) 

= - — % — en (55) 

Note that Sq(uj) decays as 1/oj at large frequencies, whereas Eq. (54) for 
classical telegraph noise decays as 1/uj 2 . 

Given these two limiting forms, the question now becomes one of how 
the SET interpolates between them. One might expect that the two results 
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should simply be added in quadrature, but as with combining thermal and 
quantum noise (see Sec. 6), this is approach is too simple. A completely 
quantum mechanical way of calculating the noise for any frequency is 
needed. This was recently accomplished by Johansson et. al [14] using an 
extension of a technique developed by Scholler and Schon [15]. Here, we re- 
derive their results using the coupled system approach outlined above. This 
method is physically motivated and allows for a heuristic interpretation of 
the final result. 

8. SET Coupled to a Qubit 

We now consider a system where the SET is coupled to a two-level system 
(i.e. a qubit), with a coupling Hamiltonian which can induce transitions 
in the TLS. Using spin operators to describe the qubit, and assuming 
operation at the degeneracy point for simplicity, where the transitions are 
fastest 10 , we have: 

H = H S et - 2^°"^ + Ao z n, (56) 

where Q is the qubit splitting frequency 11 , and A is the coupling strength. 
We can define the rates T-f and 1^ which are, respectively, the rate at which 
the qubit is excited by the SET, and the rate at which the qubit is relaxed 
by the SET. For a weak coupling (A — »■ 0), one has (c.f. Eq. 14,13 in Sec. 2): 

r l/T = ^S Q (±n) (57) 

Eq. (57) tells us that if we know the rates Tf and for a weakly coupled 
system at an arbitrary splitting frequency f2, we know the quantum noise 
Sq(Q) at all frequencies. This is the essence of the technique previously 
described, in which a qubit acts as a quantum spectrum analyzer of noise. 
Here, we mimic this approach theoretically by obtaining T| and from a 
direct analysis of the coupled system in the limit of weak coupling (A — > 0). 
The object of interest is the reduced density matrix p which describes both 
the charge n of the transistor island and the state of the qubit. We are 
interested in two quantities. First, what is the stationary state of the qubit? 
The stationary populations of the two qubit states (which are determined 

10 This amounts to maximizing the "destruction" due to the SET's noise, and the 
case where 9 — n/2, the qubit eigenstates are in the o' z — a x direction, and the SET's 
perturbation is in the —<j' x = a z direction (c.f. Eq. 42). After the noise of the SET is 
found, we can then recalculate the effects on the qubit at various n g or values of 9 by 
including the modified matrix elements in the coupling coefficient, A. 

11 Henceforth we use Q for the transition frequency, instead of the previous notation 
woi = Eoi/h, for compactness. 
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from the time-independent solution for p) will tell us the polarization of 
the qubit, and the amount of asymmetry in the noise (c.f. Eq.40). Second, 
how quickly do the qubit populations relax to their stationary value? This 
relaxation will be described by a time-dependent solution of p characterized 
by a mixing rate T m i x which is the sum of and (c.f. Eq. 41). From 
these two results we can solve for the individual values of 

To deal with the dynamics of p, we make use of the fact that sequential 
tunneling processes are completely described by lowest-order perturbation 
theory in the tunneling Hamiltonian Hp- Keeping only second order terms 
(there are no non- vanishing first order terms), one obtains the following 
standard evolution equation in the interaction picture: 



The angular brackets denote the trace over the single-particle degrees of 
freedom in the SET leads and island; pp is the equilibrium density matrix 
corresponding to the state of these degrees of freedom in the absence of 
tunneling. 12 Note that a similar density matrix analysis of a qubit coupled 
to a SET was recently discussed by Makhlin et. al [16]; unlike the present 
case, these authors restricted attention to a vanishingly small splitting 
frequency 0. 

To make progress with Eq. (58), we make a Markov approximation, 
which involves replacing pit') on the right-hand side with pit). This is 
permissible as we are interested in the slow dynamics of p. We want to find 
both the stationary solution of p, for which the Markov approximation is 
exact, and the mixing mode, a mode whose time dependence is oc e~( r f +r i)*. 
This mode is also arbitrarily slow in the weak coupling limit A — > of 
interest. Finding the stationary mode and the mixing mode correspond to 
evaluating the polarization and Tl of the qubit (c.f. Eq. 40 and Eq. 41), as 
was shown earlier for the master equation of the probabilities in Section 5. 
Note that the Markov approximation should be made in the Schrddinger 
picture, as it is in the Schrodinger picture that p will be nearly stationary 
(i.e. all oscillations associated with the qubit splitting frequency Q will be 
damped out in the long-time limit). 

Evaluation of Eq. (58) in the Markov approximation results in the ap- 
pearance of rates which are generalizations of those given in Eq. (49). Now, 
however, these rates depend on the initial and final state of the qubit- 
tunneling transitions can simultaneously change both the charge state of 
the SET island and the state of the qubit. The resulting equation is most 

12 In the diagrammatic language of Ref. [15], Eq. (58) is equivalent to keeping all (Ht) 2 
terms in the self-energy of the Keldysh propagator governing the evolution of p. 




(58) 
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easily presented if we write the reduced density matrix p in the basis of 
eigenstates at zero tunneling. For each value of island charge n, there is a 
different qubit Hamiltonian, and correspondingly a different a qubit ground 
state \g n ) and excited state |e n ). When a tunneling event occurs in the SET, 
there is a sudden change in the qubit Hamiltonian. As the qubit ground 
and excited states at different values of n are not orthogonal, tunneling 
transitions in the SET are able to cause "shake-up" transitions in the qubit. 
In the limit A — ► 0, the relevant matrix overlaps are given by: 

(9m\9n) = l~\ ( A(m fl ~ W) ) 2 = ( e m\e n ) (59) 

(e m \9n) = (60) 
Defining the frequency dependent rate r n ±i in (u;) as: 

£ 7([AiC ± i,n + M, (61) 

a=L,R 

where AE and j(AE) are defined in Eqs. (51) and (50), the required tunnel 
rates take the form: 

Em,™ = E mn (0) ~^m,n = ^m^li^) (62) 

The r + rates correspond to tunneling events where the qubit is simul- 
taneously relaxed, and thus there is an additional energy O available for 
tunneling. For large f2, tunneling processes which are normally energet- 
ically forbidden can occur if they are accompanied by qubit relaxation. 
Similarly, the T~ rates describe tunnelling events where the qubit is simul- 
taneously excited, with the consequence that there is less energy available 
for tunneling. 

Returning to the evolution equation Eq. (58), note that we do not 
need to track elements of p which are off-diagonal in the island charge 
index n- there is no coherence between different charge states, as tunneling 
events necessarily create an electron-hole excitation. Further, if we focus 
on small qubit frequencies, we may continue to restrict attention to only 
n = and n = 1 (i.e. £1 is not large enough to "turn on" tunneling 
processes which are normally energetically forbidden). Thus, there are 8 
relevant matrix elements of p- for each of the four qubit density matrix 
elements (i.e. gg, ee, ge, eg), there are two possible island charge states. 
We combine these elements into a vector p = (p gg , p ee , p gei p eg ), where 

P99 = ((°)5o|p|0,5o), (l,Si|p|l,5i)), etc. Organizing the resulting evolution 
equation in powers of the coupling A, we obtain in the Schrodinger picture: 

±p=(A + ^A 1 + §A 2 + ...)p (63) 
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We discuss the significance of the matrices Aj in what follows. 

The 8x8 matrix Ao describes the evolution of the system at zero 
coupling: 



(M 

M 

+iQ + M' 

\ o o o 



\ 



(64) 



-iQ + M' j 

with the 2x2 matrices M and M' being defined by: 



M 



-rio 

Tio 



r i 
-r i 



m' = 



-(T 

r 



10 + ^io) 



r+ + r 



10 



+ r 



10 



01 

-(r + i 



01 

+ r, 



01; 



(65) 



At zero coupling there are no transitions between different qubit states, and 
hence Ao has a block-diagonal form. There are two independent stationary 
solutions of Eq. (63) at A = (i.e. two zero eigenvectors of Ao), which 
correspond to being either in the qubit ground or qubit excited state: 



% = (Po,Pi,0,0,0,0,0,0), 



z e = (0,0,po,pi,0, 0,0,0). 



(66) 



(po,Pi) are the stationary probabilities of being in the n 
charge states: 

/ \ ( r i r 10 

(po,Pi) -- 



or n = 1 



rni + Tin Tni + r 



10 



(67) 



The existence of two zero-modes is directly related to the fact that at zero 
coupling (A = 0), the probabilities to be in the qubit ground and excited 
state are individually conserved. 

At non-zero coupling, the matrices Ai and A2 appearing in Eq. (63) 
generate transitions between different qubit states. The matrix A2 directly 
couples p gg and p ee , while Ai couples p gg and p ee to the off-diagonal blocks 
p ge and p e g. The effect of these matrices will be to break the degeneracy of 
the two zero modes of Eq. (63) existing at A = 0. After this degeneracy is 
broken, there will still be one zero mode po, describing the stationary state 
of the coupled system (the existence of a stationary solution is guaranteed 
by the conservation of probability). For weak coupling, the qubit density 
matrix obtained from po will be diagonal in the basis {| <?(„}), |e^ n ))}, which 
corresponds to the average SET charge (n) = p\. The ratio of the occu- 
pancies of these two qubit states will yield the ratio between the relaxation 
rate Tj_ and the excitation rate T|. In addition, there will also be a slow, 
time-dependent mode of Eq. (63) arising from breaking the degeneracy of 
the two A = zero modes. This time-dependent mode will describe how 
a linear combination of z g and z e relaxes to the true stationary state, and 
will have an eigenvalue A = — 1\ — IV i.e. the mixing rate. 
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Thus, we need to do degenerate second order perturbation theory in 
the coupling A to obtain the relaxation and excitation rates 1^ and P|-. 
The only subtlety here is that the matrix M is not Hermitian, implying 
that it has distinct right and left eigenvectors. Letting z represent the left 
eigenvector of Ao corresponding to the right eigenvector z, we define the 
projector matrix P as: 

P = \zg){z g \ + \z e )(z e \, (68) 

and let Pj_ denote 1 — P. As usual, degenerate second order perturba- 
tion theory requires diagonalizing the perturbation in the space of the 
degenerate eigenvectors. We are thus led to look at the matrix Q, defined 
as: 

A 2 / \ 
Q = ^2 (PA 2 P + PAiP ± [-A r 1 P±AiPj (69) 

^From the definition of Q, we may immediately identify the rates and 

T T = (~z e \Q\z g ) F t = (~z g \Q\z e ) (70) 

We thus see how the rates Tj^ arise in the present approach- they 
are related to breaking the degeneracy between two zero-modes (stationary 
solutions) which exist at zero coupling. Note that there are two distinct 
contributions to Pf^, coming from the two terms in the matrix Q: a "direct" 
contribution involving A 2 and an "interference" contribution involving Ai 
acting twice. These two terms have a different physical interpretation, as 
will become clear. 

Let us first consider the rate Pf, which describes how noise in the SET 
causes ground to excited state transitions in the qubit. For this rate, our 
approximation of only keeping two charge states will be valid for all splitting 
frequencies fl. To evaluate the "direct" contribution to this rate, which 
involves the first term in the matrix Q, note that the relevant part of A 2 
has the expected form: 

"(rs, <71) 

i.e. it consists of tunnel rates which correspond to having given up an energy 
to the qubit. Using Eqs. (69) and (70), we find: 

r T U, a - (ff - (i) 2 ( r °'^ ; .r.-, ) m 

The direct contribution to Tf has a very simple form: for each charge state 
n = 0, 1, add the rate to tunnel out of n while exciting the qubit, weighted 
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by both the probability to be in state n, and the overlap between ground 
and excited states (i.e. (A/fl) 2 ). This is very similar to how one typically 
calculates the current for a SET: one adds up the current associated with 
each charge state (i.e. a difference of rates), weighted by the occupancy 
of the state. The direct contribution to neglects any possible coherence 
between successive excitation events; as a result, it fails to recover the 
classical expression of Eq. (54) in the small-Sl limit. 

We now consider the "interference" contribution to Tf coming from 
the second term in the expression for matrix Q (c.f. Eq. (69)). After some 
algebra, we obtain the following for the interference contribution to Ty. 

2A 2 (T ) 2 

rtU = - (pdITo +Pir -i) ^ + ^ )2 (73) 

where: 

T s = r io + r io ± r oi + r oi (74) 

This contribution is purely negative, and is only significant (relative to 
the direct contribution) at low frequencies 0, < T. We can interpret this 
equation as describing the interference between two consecutive excitation 
events. For example, consider the first term in Eq. (73). This describes 
a process where a SET initially in the charge state n = undergoes a 
tunnel event to the n = 1 state, creating a superposition of qubit ground 
and excited states. At some later time the SET relaxes to the station- 
ary distribution (po,Pi) of the charge states, again partially exciting the 
qubit. Letting At represent the time between these two events, we have the 
approximate sequence: 
r~ 

\0,9o) |l,5i) + a|l,ei) 

e< nA t /2| lj ^ i)+e -inAt/2 Q | ljei) (75) 

(e inAt / 2 (3 - e- inAt / 2 a) |0, e > + ... (76) 

Here, a is the amplitude associated with qubit excitation having occurred 
during the first (n = — > 1) tunnel event, while (3 is the amplitude 
associated with excitation occurring during the second (n = 1 — > 0) tun- 
neling. These amplitudes will be given by the corresponding matrix overlap 
elements: 

A A 
a={e 1 \g )~- (3 = {e \ gi ) ~ -- (77) 

In the final state after the two tunnelings (Eq. (76)), there are two terms 
in the amplitude of the state |0, eo), corresponding to the fact that qubit 
excitation could have occurred in either the first or the second tunnel event. 
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To get a rate for this double excitation event, we should take the modulus 
squared of the final |0, eo) state amplitude, then multiply by the occupancy 
of the initial state (po) and the rate of the first tunnel event (Tio). The 
interference term in the resulting expression takes the form: 

2A 2 

T T U = (porio) x 2Re (a*Pe inAt ) = -( Po F 10 )-^- cos(ftAt) (78) 

The above expression is a function of the time At between the first and sec- 
ond tunnel events. This time is determined by the fact that the intermediate 
superposition state (Eq. (75)) corresponds to a non-stationary distribution 
of charge on the SET island, and will decay via tunneling to the stationary 
distribution (po,Pi) at a rate T^. Taking this decay to be Poissonian, and 
averaging over At, we obtain: 

T T U - -{p T 10 ) Wn2 + {T ^ 2 (79) 

This is precisely the first term in Eq. (73); the second term can be ob- 
tained in the same way, by now considering a situation where the SET is 
initially in the n = 1 charge state. As claimed, r\|j n t corresponds to the 
interference between two consecutive excitation events. The negative sign 
of this contribution can be directly traced to the matrix overlap elements 
(c.f. Eq. 77). Also, we see that the suppression of the interference term at 
large 0, results from phase randomization occurring during the delay time 
between the two excitation events. 

Returning to the total noise, we combine Eq. (73) with the direct con- 
tribution Eq. (72) to T|; comparing against Eq. (57), we obtain the final 
expression for Sq(Q) at all negative frequencies: 

s Q (-\n\) = PoiYo + Piroi (80) 

^ 2 + i(rro + r5 ) + r - 1 + r + 1 ) 2 

Note for large \Q\ (i.e. > max(A£7gi, A£f ) ~ V DS /2), Sq(-Q) will 
vanish identically at zero temperature. Physically, this cutoff corresponds 
to the largest amount of energy the SET can give up to the qubit during a 
single tunnel event; giving up more energy would suppress the event com- 
pletely (i.e. the tunnel rates have a step-function form at zero temperature, 
c.f. Eq. (50)). If one included higher order processes in the tunneling (i.e. 
went beyond sequential tunneling), correlated tunneling events involving 
the full voltage drop over both junctions, Vds, would move this cutoff to 
higher values of absolute frequency. 

We now turn to the relaxation rate r\ , and hence the positive frequency 
parts of Sq. The calculation proceeds exactly as that for F^, the only 
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Figure 5. Symmetrized SET charge noise as a function of frequency, for typical SET 
parameters: 3 = 1, E c /k B = 2K, TV = 0.33, eVbs = E c , and T = 2Qm.K. The dashed 
line is the classical telegraph noise (Eq. (54)), while the dot-dashed line is the noise of 
two parallel tunnel junctions (Eq. (55)). Inset: full (non-symmetrized) quantum noise for 
identical SET parameters; the dashed line is the symmetric classical telegraph noise. 



modification being that one now needs to include the charge states n = 2 
and n = — 1, as the SET could absorb enough energy from the qubit to 
make transitions to these states possible. We can combine the result for 1^ 
with Eq. (80) to obtain a single, compact expression for the noise at all 
frequencies first obtained by Johansson et. al [14]: 13 

s / x = pq FioH + r-i,oM] +pi [ToiM + r 2 iM] . . 
^ + i[r 1 oH + r 1 o(-^) + roiM + r i(-^)] 2 

Shown in Figure 5 is the symmetrized noise Sq(lo) + Sq(-lo) for typical 
SET parameters. One can clearly see abrupt changes in the slope of this 
curve; each of these kinks corresponds to a threshold frequency at which a 
given tunneling process either turns on or turns off. For comparison, curves 

13 Eq. (81) ignores additional order g/(2ir) terms which arise in the denominator at 
positive frequencies large enough to turn on tunneling to higher charge states; such terms 
are clearly negligible in the sequential tunneling regime due to the smallness of g/(2n). 



rjsdraftl2.tex; 1/02/2008; 22:49; p. 28 



29 




gate 1 'gate 



Figure 6. a) Average charge state of a Cooper pair box coupled to a SET, as a function 
of box gate voltage, using identical SET parameters as above. The box parameters 
are Ec/kB = 0.5K, Ej/Jcb = 0.25K and the coupling constant is k = 0.04. We also 
include relaxation effects due to a 10% coupling to a 50f2 environment. The dashed curve 
corresponds to assuming the SET produces classical telegraph noise, the solid curve 
corresponds to using the full quantum noise of the SET, and the dashed-dot curve is the 
box ground state, b) The relaxation time T\ for the same system, as a function of box 
gate voltage. 



corresponding to classical telegraph noise and to the uncorrected noise of 
two tunnel junctions are also shown. At low frequencies the symmetrized 
true noise matches the classical curve; for higher frequencies, it lies above 
the classical curve but below the curve corresponding to the uncorrelated 
case. The inset of this figure shows the both the negative and positive 
frequency parts of Sq(uj). 

It is easy to check that in the limit u — ► 0, Eq. (81) recovers the classical 
telegraph expression of Eq. (54). In the high-frequency, zero temperature 
limit, one can also see that Eq. (81) approaches the uncorrelated result of 
Eq. (55) from below: 

s Q {u>) en \ 2 ; e( w )-£ (82) 

Note that at high frequencies, it is only the "direct" terms which contribute 
to the noise- the interference contribution is not important in the limit 
of uncorrelated tunneling. The fact that the noise approaches the high 
frequency limit from below results from the tendency of charging energy 
induced correlations (which are present for a finite uj/Eq) to suppress 
fluctuations of n, and thus suppress the noise. Note that the interpolation 
between the low and high frequency limits here is very different than, e.g., 
interpolating between thermal noise and zero-point fluctuation noise in a 
tunnel junction. In the latter case, one is effectively combining two sources 
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of noise; here, one is simply turning off correlations brought on by the 
charging energy by increasing u. 

Finally, shown in Figure 6 is the average charge state of a Cooper pair 
box coupled to a SET with identical parameters to that in Fig. 5. We have 
also included here the relaxation effects of the environment, modelled as 
in Sec. 6 as a 50f2 impedance. Note that even near the box degeneracy 
point, there are large deviations between the result obtained using the full 
quantum noise of the SET and that obtained from using only classical 
telegraph noise. In Fig. 6b, we show the relaxation rate T\ for the same 
system. Note that the differences between using the full quantum noise and 
the classical expression are not so evident here. 

9. Summary 

In this article, we have emphasized the need to discuss quantum noise 
processes using their two-sided spectral densities. Because of the quantum 
nature of noise, the positive and negative frequencies are generally unequal, 
in order to account for spontaneous emission. A two-level system was shown 
to be an ideal spectrum analyzer for probing the quantum nature of a 
noise process or reservoir. With the advent of real electrical circuits which 
behave as coherent two- level systems (e.g., [5], [7]), we can now build and 
use quantum electrical spectrum analyzers. We also described the use of 
a qubit as a theoretical tool, by following the evolution of the density 
matrix of a TLS coupled to the noise-producing system of interest. This 
technique appears to be quite powerful, as it can yield analytical results 
for the full quantum noise spectrum of a wide variety of devices, including 
the superconducting SET [17]. The distinction between the classical noise 
and the quantum noise, found in this way, leads to dramatically different 
predictions (c.f. Fig. 6 and Ref. [17]) for continuous measurements of qubits 
with an SET. The "coupled-system" calculational approach also allows 
predictions of the dephasing by the measurement, the performance relative 
to the Heisenberg uncertainty limit [13], the fidelity of single-shot measure- 
ments of the qubit states, and the effects of strong coupling to the qubit. 
The combined theoretical and experimental advances raise many interesting 
possibilities for testing our understanding of quantum measurement theory 
with mesoscopic devices. 
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